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Abstract 


We calculate expected event rates for direct detection of relic neutralinos 
as a function of parameter space of the minimal supergravity model. Numer¬ 
ical results are presented for the specific case of a 73 Ge detector. We find 
significant detection rates (R > 0.01 events/kg/day) in regions of parame¬ 
ter space most favored by constraints from B —> X and the cosmological 
relic density of neutralinos. The detection rates are especially large in re¬ 
gions of large tan /3, where many conventional signals for supersymmetry at 
collider experiments are difficult to detect. If the parameter tan /3 is large, 
then there is a significant probability that the first direct evidence for super- 
symmetry could come from direct detection experiments, rather than from 
collider searches for sparticles. 
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I. INTRODUCTION 


Luminous matter comprises less than = pjp c ~ 0.01 of the matter density of 

the universe [|I]|§, where p c = is the critical closure density of the universe, and 

H = lOOh km/sec/Mpc is the scaled Hubble constant, with 0.5 < h < 0.8. In contrast, 
most inflationary cosmological models require p = p c , i.e. a flat universe. The difference 
in these matter density values can be reconciled by hypothesizing the existence of dark 
(non-shining) matter (DM) in the universe. Experimental evidence for galactic dark matter 
comes from the enclosed mass versus velocity plots measured for clouds of neutral hydrogen 
rotating about galactic centers (galactic rotation curves), which imply D > 0.03 — 0.1. An 
understanding of galactic clustering and galactic flows points towards even larger values 
of ff ~ 0.2 — 1, possibly in accord with inflationary cosmology. Calculations of Big Bang 
nucleosynthesis can only allow the baryonic contribution to the matter density of the universe 
to be ff baryonic ~ 0.01 — 0.1, so that if p = p c , the bulk of dark matter in the universe must be 
non-baryonic. Some candidates for non-baryonic dark matter from particle physics include 
neutrinos with eV scale masses (hot dark matter: HDM) and WIMPs (weakly interacting 
massive particles), such as axions or the lightest neutralino in supersymmetric (SUSY) 
models. Some compelling models of structure formation in the universe which take into 
account COBE measurements of the anisotropy in the cosmic microwave background actually 
prefer a “mixed dark matter” (MDM) universe, with ~ 10% baryonic matter, ~ 30% HDM, 
and ~ 60% cold dark matter (CDM, comprised of WIMPs). 

In this paper, we focus attention on the lightest SUSY particle, or LSP— usually the 
lightest neutralino— from supersymmetric models. At very early times after the Big Bang, 
neutralinos would have existed in thermal equilibrium with the primordial cosmic soup of 
particles and radiation. As the universe expanded and cooled, temperatures dropped so 
low that neutralinos could no longer be produced, although they could still annihilate away. 
Ultimately, the expansion rate of the universe exceeded the annihilation rate so that a 
relic density of neutralinos would be locked in. Numerous estimates of the neutralino relic 
density Qh 2 as a function of SUSY model parameter space have been made ||f|]. Models that 
predict ilh 2 < 0.025 cannot even account for the dark matter needed for galactic rotation 
curves, while values of f lh 2 > 1 would yield a universe younger than 10 billion years old, in 
contradiction at least with the ages of the oldest stars found in globular clusters. 

A consequence of the SUSY dark matter hypothesis is that a non-relativistic gas of 
neutralinos fills all space. To test this hypothesis, a number of direct detection experiments 
have been built or are under construction |J. The general idea behind these experiments 
is that relic neutralinos (or other possible WIMPs) could scatter off the nuclei in some 
material, depositing typically tens of keV of energy. Some examples of how the thermal 
energy could be detected include: i.) via changes in resistance due to a slight temperature 
increase (bolometry), ii.) via a magnetic flux change due to a superconducting granule phase 
transition, or in.) via ionization. The technical challenge is to build detectors that could 
pick out the relatively rare, low energy neutralino scattering events from backgrounds mainly 
due to cosmic rays and radioactivity in surrounding matter. Future detectors are aiming to 
reach a sensitivity of 0.1-0.01 events/kg/day. In this way, the first evidence for SUSY might 
come from direct neutralino detection rather than from accelerator experiments. 

Since the pioneering paper by Goodman and Witten j6j, calculations of neutralino- 


2 



nucleus scattering have seen continual improvements The first step involved in a 

neutralino-nucleus scattering calculation is to calculate the effective neutralino-quark and 
neutralino-gluon interactions. The neutralino-quark axial vector interaction leads in the 
non-relativistic limit to a neutralino-nucleon spin-spin interaction, which involves the mea¬ 
sured quark spin content of the nucleon. To obtain the neutralino-nucleus scattering cross 
section, a convolution with nuclear spin form factors must be made. The neutralino-quark 
and neutralino-gluon interactions (via loop diagrams) can also resolve into scalar and tensor 
components. These interactions can be converted into an effective scalar neutralino-nucleon 
interaction involving quark and gluon parton distribution functions. A neutralino-nucleus 
scattering cross section can be obtained by convoluting with suitable scalar nuclear form fac¬ 
tors. The final neutralino detection rate is obtained by multiplying by the local neutralino 
density (estimates are obtained from galaxy formation modeling), and appropriate functions 
involving the velocity distribution of relic neutralinos and the earth’s velocity around the 
galactic center. 

In this paper, we calculate event rates for direct detection of relic neutralinos left over 
from the Big Bang. For illustration, we present detailed calculations for a ' 3 Ge detector; 
a 73 Ge detector has a sizable nuclear spin content J — |, so that it would be sensitive to 
both spin and scalar neutralino-nucleus interactions JD). Our goal is to compare the reach 
for supersymmetry by direct detection experiments with constraints on model parameter 
space from relic density calculations, collider search limits, and limits from B —»■ X s j mea¬ 
surements. In addition, we will compare the reach of direct detection experiments with the 
previously calculated reach of collider facilities such as LEP2, the Fcrmilab Tevatron Main 
Injector (MI) upgrade and the CERN LHC pp collider. 

We work within the framework of the paradigm minimal supergravity (mSUGRA) model 
nj. This model assumes the minimal supersymmetric standard model, or MSSM, is valid 
at all energy scales from M wea k up to Mqut — 2 x 10 16 GeV. The mSUGRA model could 
arise as the low energy limit of a supergravity theory, where supersymmetry is broken in 
the hidden sector of the model at energy scale M ~ 10 10 GeV. Supersymmetry breaking is 
communicated to the observable sector via gravitational interactions, leading to soft SUSY 
breaking mass terms of order the electroweak scale, m ~ 100 — 1000 GeV. At the GUT scale 
(with the added assumption of an approximate global U{n ) symmetry for the mSUGRA La- 
grangian), this leads to a common mass for all scalars mo and a common trilincar coupling 
A 0 . Motivated by the apparent unification of gauge coupling constants, it is also assumed 
that all gaugino masses are unified to mi/ 2 at M GUT . The weak scale sparticle spectrum is 
derived from renormalization group (RG) running of the SUSY soft breaking parameters. 
Requiring radiative electroweak symmetry breaking allows the determination of the super¬ 
potential Higgsino mass squared p 2 , and allows the expression of the soft SUSY breaking 
bilinear term B in terms of tan (3 , the ratio of vev’s of the two Higgs fields. Thus, all sparticle 
masses and couplings are derived in terms of the parameter set 


mo, mi/ 2 , A 0 , tan/i, and sign(p). 

We adopt the mSUGRA spectrum calculation encoded into the event generator IS A JET 

7.29 |lj . 


In Sec. 2, we present details of our calculation for neutralino scattering events off a 73 Ge 
detector. In Sec. 3 we present numerical results for event rates in the parameter space of the 


3 





mSUGRA model, and compare to neutralino relic density contours. In Sec. 4, we compare 
our results to B —> Xg'y constraints, and to collider search constraints. Finally, we compare 
the search capability of a 73 Ge dark matter detector to the reach of various future collider 
experiments for mSUGRA. We find usually that a DM detector attaining a sensitivity of 
0.01 events/kg/day will have a greater reach into mSUGRA parameter space than either the 
LEP2 or Tevatron upgrades, via their searches for sparticles. However, such dark matter 
detectors can only probe a fraction of the parameter space that gives rise to a reasonable 
relic density. A complete exploration of the cosmologically interesting mSUGRA parameter 
space will have to await the CERN LHC pp collider, or an e + e^ or g + g~ collider operating 
at y/s ~ 1 TeV. In Sec. 5, we give a summary and some conclusions. 


II. CALCULATIONAL DETAILS 


A. Dark Matter Detection: Theory 


The effective elastic scattering Lagrangian can generally be divided into two parts: 


ye// _ ye// I ye// 
^elastic ^scalar ' spin• 


( 2 . 1 ) 


We examine first the scalar Lagrangian, which receives contributions from neutralino-quark 
interaction via squarks and Higgs bosons exchange, and from neutralino-gluon interactions 
at one-loop level involving quarks, squarks and Higgs bosons in the loop diagrams. On the 
parton level it is expressed at scale Q (typically ~ ) as [§] 


Cellar = fqZlZiqq + 9q 


-2iZi'y fl d u Z 1 Q i2)flu - -rriqrn^ZiZiqq 
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Here 
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are traceless twist-2 quark and gluon operators, and for the sake of brevity the effective 
couplings are given in the Appendix. Using nucleonic matrix elements 
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(here, q 2 is the momentum transfer squared, and is the nucleon four-momentum) and 
introducing the parton distribution functions qjy, q jv and g at, and using 


(N\F^F^(q 2 - 0)|Al) = 


4 









(N\m g qq\N) = m N ff q \ 


it is possible to convert Lagrangian (2.15) into an effective neutralino-nucleon Lagrangian 


Cellar = fpZlZ^^p + fnZ&Vrf&r 
Evaluation of /jv for N = p,n yields 
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where the various coupling constants are given in the Appendix. Here we have used the 
general definition of an n-th integral moment f(n, Q ) = Jq x n ~ 1 f(x , Q) dx and applied it to 
the parton distribution functions. The differential cross section for a neutralino scattering 
off a nucleus X% with mass is then expressed as 


da' 


scalar 


d\q\‘ 


7 TV 


:[Zf p + (A-Z)f n YF\Q r ), 


(2.5) 


where q = 


m A m~ 

£] 


-v is the transfered momentum, Q r = and F 2 (Q r ) is the scalar nuclear 


form factor. 

Interaction between the neutralino and quark spins is described by a parton-level La¬ 
grangian 0 


£ e s f J n = dqZiYlsZ^^q 

which translates with the help of nucleonic spin matrix elements 

< A , |? 7 p 75 ?|iV >= 2 s„A </ (a ' 1 


( 2 , 6 ) 


into 

C e /Jn = 2 V2 {(ipZ^^Zi^pS^Vp + a n Zi7 /i 7 5 Z 1 T n s At T n ) , (2.7) 

explicitly involving the nucleon spin vectors s M . Coefficients 

«p = 4E d q Aq (p \ a n = d q Aq {n) (2.8) 

V 2 u ^ s V 2 u d s 

depend on experimental values of A q ^, which are affected by significant uncertainties which 
lead to variations in the cross section. More details on the couplings are again given in the 
Appendix. For a nucleus with total angular momentum J, the spin interaction differential 
cross section takes the form 
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d a spin 

d\q\ 2 


^A 2 J(J + 1) 
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sm 

S{ o) ’ 


(2.9) 


where is the nuclear spin form factor normalized to 1 for pointlike particles, and 

A = j| a p (S p ) +a n (S n )]. The quantities ( S p ) and (S n ) represent the expectation value of the 
proton (neutron) group spin content in the nucleus. 

Putting both scalar and spin interaction contributions together and convoluting them 
with the local neutralino flux (which depends on the local relic density p% ), the differential 
detection rate is calculated to be 

= -^=-^-T(Q r )\[Zf r +(A- Z)f n ] 2 F 2 (Q T ) 

dQr V 7T 3 m^V 0 l 

+8A 2 J(J + 1)^M}, (2.10) 


where r>o ~ 220km/s is the circular speed of the Sun around the center of our galaxy and 

T(Qr) = ^ [°° fj ±dv ( 2 . 11 ) 

" J v min U 

integrates over the neutralino velocity distribution. 


B. Detection Rates in Germanium Detectors 


In order to obtain a quantitative estimate of neutralino detection rates in the mSUGRA 
framework which could be compared to phenomenological constraints and collider reaches, 
we have evaluated the rate for the case of a ' 3 Ge detector. In our calculations, we have taken 
the neutralino density to be the local relic density consistent with galactic formation models 
including varying baryonic mass fraction in the galaxy |L2|. In the simplest Gaussian model 
taking into account the motion of the Sun and Earth, the integrated velocity distribution 
can be written as [01 


T(Qr) = 


Vkvq 

4m 


Erf( Vmin + Ve ) - Erf( Vmin Ve ) 
Vo V 0 


( 2 . 12 ) 


lQr(m~ +m . A ) 2 

with v rnin = ./——, and where the Earth velocity v e is given by 
V m Zi m A 


Ve = Vo 


1.05+ 0.07 cos 


,27T (t — t p ) 
1 yr 


(2.13) 


with t p ~ June 2. The general nuclear properties needed are the 73 Ge mass, rric e = 67.93 
GeV, and its total spin J — |. 

To determine the scalar contribution, it is necessary to compute the parton distribution 
integrals at a scale defined by the average squark mass and neutralino mass Q 2 ~ ( m~—m 


zd 


We employ the CTEQ3L parton distribution function parametrization JT3[] for numerical 
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calculation. The most recent values of the matrix element coefficients 
1 — fr!) were compiled by 


AN) 
J Tq 


and f$) = 


J Tq 

Ap) 


giving #1 = 0.019, ,/?j = 0.041, = 0.023, = 0.034 


and flpj = fxs = 0-14. We adopt the Saxon-Woods scalar form factor suggested in |L4 


F(Q r ) = 3jl ( qRo K -fa*)\ 
qR\ 

where Ri = \JR 2 — 5s 2 , R— A 3 x 1.2 fm, ji is a spherical Bessel function and s = 1 fm. 

In the case of the spin interaction, the spin analogues of the parton distibution functions 
are much less well known, and we take (T5|] A u p = A d n = 0.78, A d p = A u n = —0.5 and 


A s p = A s n = —0.16. Theoretical predictions for the spin content of the two nucleon groups 
in the nucleus and the spin form factor are very model dependent. To make a consistent 
choice, we follow [|16j], where for ' 3 Ge ( S p ) = 0.03, ( S n ) = 0.378 and the form factor is given 
as 


S(q) — ( a p + a n ) 2 Soo(q ) + ( a p — a n ) 2 Sn(q ) + ( a 2 — a^)Soi(q) 

and the Sij individual form factors are evaluated as polynomial fits to data. 

As a final step in the calculation, the differential rate -j§- must be integrated over Q r 

ClLJr 

ranging typically from ~ 0 to less than 100 keV. We will show the rate R in events/kg/day 
for a particular choice of the assumed local neutralino density pg = 5 x 10 -25 gcm -3 ||12|| . 
The standard lore is that the local halo density uncertainty should be roughly a factor of two; 
a more accurate prediction based on galaxy formation models would be needed to reduce it. 


III. NEUTRALINO DETECTION RATES IN MSUGRA PARAMETER SPACE 


We show our first numerical results for direct detection of neutralinos in Fig. 1, where 
we plot contours of scattering events/kg/day for a 73 Ge detector as a function of mSUGRA 
parameters m 0 vs. rn 1 / 2 with A 0 = 0, tan/? = 2 and a) /i < 0 and b) p > 0. The region 
labelled TH is excluded by theoretical consideration: either the LSP is charged or colored 
(not the lightest neutralino), or radiative electroweak symmetry breaking is not properly 
attained. The region labelled EX is excluded by collider searches for SUSY particles. By 
far the strongest of these for mSUGRA is the recent limit from LEP2 that > 85 GeV 


a We also show the contour of constant Qh 2 = 1 (solid); beyond this contour, mSUGRA 


parameter space points lead to universe with age less than 10 billion years. In addition, the 
dot-dashed contours correspond to Qh 2 = 0.15 and 0.4; the region between these contours 
is favored by MDM cosmological models. 

We find in frame a) that the direct dark matter detection rates are uniformly low through¬ 
out the parameter space shown; for all points sampled, we found R < 0.01/kg/day, which 
is less than the goal for such detectors at least in the near-term future. In frame b), how¬ 
ever, for p > 0, we find larger rates for dark matter detection, with a significant fraction 
of parameter space with m 1/2 < 200 GeV accessible to dark matter detectors achieving a 
sensitivity of R > 0.01/kg/day. Note that the R > 0.01/kg/day region overlaps with the 
lower portion of the region favored by a MDM universe. 

I 11 Fig. 2, we show similar results, except now we take tan/3 = 10. In Fig. 2, the region 
below the dotted contour is where Qh 2 < 0.025— too small to account for the galactic 
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rotation curves. In frame a), we find considerably larger detection rates than for Fig. la, 
with R reaching values ~ 0.1 in the lower-left. However, in this region, flh 2 < 0.025 so that 
the highest detection rates exist in an uninteresting region of relic density. DM detectors with 
sensitivity R > 0.01/kg/day can probe just the lower portion of the MDM cosmologically 
favored region. Likewise, in frame b ), event rates are larger as well than in Fig. 16, with 
the entire region shown below m 1/2 — 200 GeV accessible to DM detectors able to achieve 
a counting rate of R = 0.01/kg/day. 

Fig. 3 is similar in construction to Figs. 1 and 2, but with tan/3 = 35. Note, however, 
the expanded scale relative to Figs. 1 and 2. The first peculiarity of note is that the TH 
region has expanded considerably. This is due to the magnitude of the r Yukawa coupling at 
large tan (3, which drives the fi mass to lower values than corresponding sleptons from the 
first two generations. The expanded TH region in the upper left is thus where m^ < m-g, 
so that the f\ is the LSP instead of the lightest neutralino. 

I 11 Fig. 3 a and b, we see that the regions of cosmologically interesting relic density flh 2 
are much larger than for the low tan (3 cases. In this case, the enhanced Higgs coupling to 
bb and rr at large tan f3 gives rise to a very broad resonance structure so that s-channel 
annihilation of neutralinos is possible over a very large region of parameter space. Note in 
particular that no upper limit on SUSY particle masses is evident in these plots from the 
Qh 2 < 1 constraint. 

We find in Fig. 3 that the dark matter detection rate R has grown even larger for 
low mi /2 values, relative to Figs. 1 and 2, so that R exceeds 1 event/kg/day in the lower 
left! I 11 this region, however, Vth 2 < 0.025 so that again the largest detection rates are in 
a cosmologically uninteresting region. We see that R can exceed values of 0.01/kg/day for 
m 1/2 as high as 300 GeV, which corresponds to a reach in rrig of ~ 750 GeV! 

The large increase in DM detection rate for large tan/3 has been noted previously by 
Drees and Nojiri [||. This is a fortuitous result: the region of parameter space where DM 
detection is easiest is precisely the region of parameter space where collider detection of 
SUSY particles is most difficult (this will be discussed in more detail in Sec. 4). We show 
the explicit variation in R with tan (3 for the mSUGRA point m 0 , mi/ 2 , A 0 = 150, 200, 0 GeV 
in Fig. 4. Here we see the event rate increasing by more than one (two) orders of magnitude 
with tan/3 for /1 > 0 (/r < 0). Fig. 4 has some resemblance to Fig. 3 of the second paper 
of Ref. [|]; in our case (for a 73 Ge detector instead of a spin-0 ,6 Ge detector), the large 
spin-spin interaction causes the main difference between the plots. 

For 73 Ge, the spin coupling is large, and can cause the axial-vector interaction to exceed 
the scalar interaction even for nuclei as heavy as Ge. This is shown in Fig. 5, where we 
plot the ratio R sp i n /Rscalar for Ike same mSUGRA point as in Fig. 4. Here, we see that 
the spin interaction actually dominates for tan/3 ~ 6, for /i < 0. This is contrary to naive 
expectations that the spin interaction is always sub-dominant for nuclei with A > 20. 

In Fig. 6, we show an example of the variation in dark matter detection rate versus 
variation in the trilinear soft-breaking term A 0 , for mo, Tn\n — 150, 200 GeV, and tan f3 = 10. 
From this plot, we see that the DM detection rate can vary with A 0 by an order of magnitude; 
this sort of variation must be taken into account if dark matter detection experiments ever 
try to obtain limits on mSUGRA parameter space. 

Finally, we note that a seasonal variation in the DM detection rate is expected ||. This 
is caused by accounting for the earth’s velocity about the sun, while at the same time 
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accounting for the sun’s velocity about the galactic center. In Fig. 7, we plot the DM 
detection rate as a function of month of the year (beginning with Jan. 1). The maximum 
detection rate occurs around June 2, although the seasonal variation amounts to less than 
a per cent, so that very high counting rates would be necessary to detect this. We note as 
well that if DM detectors are sensitive to the direction of collision products, these should 
also depend on the time of day and season of the year. 


IV. COMPARISON WITH SUSY SEARCHES AT COLLIDERS 


In this section, we compare our results for direct detection of neutralino DM with con¬ 
straints from B —> X s j searches, and with expectations for various collider searches. In 
these comparisons, we are restricted to values of tan/? < 10, for which detailed collider and 
B —> Xg'y calculations are available. Only recently has the event generator ISAJET been 
upgraded to handle large tan/? cases, so detailed studies for collider expectations still need 
to be made 0- 

The rare decay B —> X s j has been shown to yield rather strong constraints on super- 
symmetric models, due to loop amplitudes containing charginos, neutralino, gluinos, squarks 
and charged Higgs bosons ]17|. In a recent paper pU| , a QCD-improved calculation of the 
B —> Xg'y branching fraction has been made for the mSUGRA model. The sensitivity to 
variations in the QCD renormalization scale has been reduced considerably compared to 
previous results. In this paper, only chargino, charged Higgs and JU-boson loops have been 
included, which is appropriate for small to moderate values of tan/?. Comparison of cal¬ 
culated branching fractions to recent results from the CLEO experiment has resulted in 
identification of regions of mSUGRA model parameter space which are excluded at 95% 
CL. For large tan/?, g and Z,- loops will also be relevant (see Borzumati, Drees and Nojiri, 
Ref. JT9[]). In this case, the derived constraints will depend sensitively upon details of the 
assumed structure of high scale squark mass matrices, so that the implications will be much 
more model dependent. 

In Fig. 9, we again show the m 0 vs. m\j 2 plane for tan/? = 2. This time, we show in 
Fig. 9a the region excluded by the 95% CL CLEO result compared to mSUGRA model 
calculations. We exclude parameter space points where the B —> Xg'y branching ratio falls 
outside the CLEO 95% CL limits for all choices of renormalization scale ^ < Q < 2 mt,. 
The excluded region is in the lower left of frame a), where the dark matter detection rate 
is largest, although still smaller than R = 0.01/kg/day. In this frame, we also show the 
region that can be searched by the LEP2 e + e _ collider operating at yfs = 190 GeV, and 
accumulating ~ 500 pb _1 of integrated luminosity (dotted curve). The left-side bulge in this 
contour is where LEP2 is sensitive to selectron searches, while below the right-hand side, 
which asymptotically approaches 2 — 100, is where LEP2 is sensitive to chargino pair 

a 


searches 


We also note that LEP2 is sensitive to almost the entire plane shown via the 


e + e _ —> Zh search channel; in this case, however, it will be difficult to tell the light SUSY 
Higgs h from the SM Higgs boson. Finally, we also show the dashed contour, which is the 
reach of the Tevatron Main Injector (MI: y/s = 2 TeV; integrated luminosity = 2 fb _1 ). This 
latter curve mainly results from regions where the clean trilcpton signal from W\Z 2 —■► 3£ is 
observable above SM backgrounds E3- By comparing all the contours of Fig. 9 a, we note 
that the low mi/ 2 region of the MDM-favored region is excluded by B —> Xg'y, which also 
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excludes much of the region open to SUSY discovery by LEP2 and Tevatron MI searches. If 
mSUGRA is correct, with tan (3 — 2 and p < 0, then LEP2 may well discover the light Higgs 
scalar h, but probably there will be no direct detection of SUSY DM and SUSY discovery 
will have to wait for the CERN LHC pp collider, which has the ability to explore the entire 
SUSY parameter space with relatively low luminosity [|23| . 

In Fig. 96, none of the mo vs. m \/2 plane is excluded by B —> X s 7 . The reach of LEP2 
is almost always below the MDM-favored region. However, both the Tevatron MI as well as 
DM detectors sensitive to R ~ 0.01/kg/day can explore the lower limits of the MDM region. 
In this case, the MI has somewhat of a better reach for SUSY than direct DM detection 
experiments, while LEP2 would stand again a good chance to find the light Higgs scalar h 

& 

I 11 Fig. 10 a, we show the plane for tan (3 = 10 and p < 0. The outstanding feature here is 
that B —> Ad. 7 data exclude almost the whole plane below m ^ 2 ~ 350 GeV (corresponding 
to < 900 GeV), including essentially all of the MDM-favored region [P0| ! If mSUGRA 
exists in this parameter plane, then direct DM detection experiments, LEP2 and Tevatron 
MI will find no evidence of SUSY, and SUSY discovery will have to await the LHC (although 
B —> Xg'y experiments may have a strong hint of new physics). 

Fig. 106 shows the mo vs. plane for tan/? = 10 and p > 0. In this case, none of 
the plane shown is excluded by B —> Ad. 7 . In fact, the mSUGRA region around — 200 
GeV gives a better match to CLEO data than does the SM! For this case, the reach of LEP2 
is below the MDM region, and in addition the Higgs scalar h is too heavy for detection 
at LEP2. The Tevatron MI can explore a portion of the MDM favored region. However, 
a DM detector with sensitivity R ~ 0.01/kg/day will be sensitive to the additional region 
with mi / 2 /$ 200 GeV, and m 0 > 150 GeV, which is inaccessible to Tevatron or LEP2 
searches. This region includes the entire MDM-favored s-channel annihilation corridor for 
which mi /2 ~ 180 GeV. This latter region is also favored by the B —> A /,7 measurement 

[23. 

Finally, we content upon the large tan f3 region, for which detailed collider studies and 
B —> Ad. 7 calculations have yet to be made. Drees and Nojiri pointed out that as tan f3 in¬ 
crease, the Higgs pseudoscalar mass decreases significantly, so that over much of param¬ 
eter space the neutralino relic density decreases, mainly due to s-channel ZZ —> A, H —> 66 
annihilation reactions |3j]. We have verified this with the relic density contours presented 
in Fig. 3. I 11 a recent paper |1|, detailed calculations of sparticlc masses, production and 
decay processes at large tan/3 were reported. It was noted that for large tan (3, the W\ 
and Z 2 branching ratios to r leptons and 6 -quarks increases due to Yukawa coupling effects, 
which leads to a diminution of the corresponding branchings to easily detectable e and p 
states. This generally ought to make SUSY detection much more difficult for the Tevatron 
MI than corresponding mSUGRA points with low tan/?. However, we note from Fig. 3 that 
the large tan/? region is precisely where DM detection rates can be largest. Hence, if tan/? 
is large, it is possible that the first evidence for SUSY might come from direct dark matter 
detection, rather than from LEP2 or MI searches. It is expected that the CERN LHC pp 
collider will still be able to cover the entire mSUGRA parameter space even if tan/? is large, 


at least via multijet+/? T searches [|18 
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V. SUMMARY AND CONCLUSIONS 


In this paper, we have presented expected event rates for direct detection of neutralinos 
by cryogenic dark matter detectors. To be specific, we chose to perform calculations for a 
,3 Ge detector. Many other choices of materials are possible, and analogous calculations can 
be made by using different nuclear A, Z and J values, and by using different nuclear form 
factors. Our main results were presented in Figs. 1-3 as functions of mSUGRA parameter 
space. When interpreting these results, the intrinsic uncertainties in the calculations should 
be kept in mind. Roughly, we expect a factor of ~ 2 uncertainty from each of i.) calculations 
in leading-log QCD, ii.) uncertainty in knowledge of the local dark matter relic density 
and in.) uncertainty in nuclear form factors and in quark contributions to nucleon spin. 
Simplistically adding these in quadature would imply our results for event rates are only 
reliable to a factor 3-4. In addition, variation of the parameter Ao can cause changes in the 
DM detection rates by up to an order of magnitude (see Fig. 6). 

One of our main goals in this paper was to compare the reach for supersymmetry by 
dark matter detectors against the reach for SUSY via collider experiments and rare decay 
searches. A dark matter detector reaching a sensitivity of R ~ 0.01/kg/day usually would 
have a better reach in mSUGRA parameter space than LEP2 would for SUSY particles 
(but not if one includes LEP2 sensitivity to Higgs bosons). Dark matter detectors can be 
comparable to the Tevatron MI in terms of reach for mSUGRA for low values of tan/3. 
However, for large tan /3, dark matter detectors have an increased event rate, whereas the 
reach of the Tevatron MI will likely diminish relative to capabilities at modest values of 
tan (3. For the large tan /3 case, the first direct evidence for SUSY may well come from direct 
detection experiments. We note that even if dark matter detectors only achieve a sensitivity 
of 0.1/kg/day, they would still have a substantial reach for SUSY in the large tan/3 region. 
The relative capabilities at low and high tan /3 underscores another facet of complementarity 
between collider search experiments, and direct detection of dark matter. Obviously, if one 
detects SUSY at collider experiments, it would still be fascinating to verify the existence 
and properties of neutralino dark matter. 
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APPENDIX A: EFFECTIVE COUPLINGS FOR NEUTRALINO-NUCLEON 

SCATTERING 

In this Appendix we summarize effective couplings entering the scalar and spin La- 
grangian for neutralino scattering on nucleons, and express them in terms of the MSSM 
coupling constants and masses. Most of the time we follow reference ||. 

Both interactions are mediated by squarks, Higgs bosons or the Z° boson connecting 
either directly or through a loop the neutralino lines of propagation with quark (gluon) lines 
in the nucleon. The lightest neutralino-quark-squark Lagrangian reads 

C z tq q = E Z t (a~ + b-jsjqq, + h.c. (A.l) 

<= 1,2 

and the a and b coefficients take the form 

“5 = ^{~mee,N n + - e,^0 w ) + £^-N ,J 

+ AfeleeX - -JL-N B e,sm>e w - f^j}, (A.2) 

h = - e,^0 w ) - ^-N„] 

+ M i2 [~ee q N' n H- ^—N' a e q sin 2 6 W - | ’—Nix] \, (A.3) 

cos t?w zmw ) 

where 

(q 1 \ = (M 11 M 12 \ ( Ql\ ( cos 9 q sin 6 q \ (q L \ ( . . 

\qj \M 21 M 2 J’WJ v — sin 9 q cos d q )\q R )’ { 

and x = 3(4) for a down (up) type quark. The 4 x 4 N matrix diagonalizes the neutralino 
mass matrix and 


( N ji \ = ( cos °w sin 9 W \ / Njt \ 
\N’ j2 ) V — sin dw cos 9 W ) \N j2 ) 

The f q coupling in (2.15) can be split into two parts 

f — f(5) 4 . f( H ) 

Jq Jq ~ Jq ) 


where the squark part is 


/® 




qt 



- (m^ + m q ) 2J 


and the Higgs exchange part is 




= m, 


c d) c (i) 
z 1 


q L—t 

3=1,2 


mV 


(A.5) 


(A. 6 ) 


(A.7) 


(A. 8 ) 
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Mixing in the Higgs sector results in 


c~ } = - (gN e2 - g'N n ){N t3 sin a + A^ 4 cos a) 
for the lighter CP-even Higgs and 

c~ = ^(s N t 2 - g'N n )(N, t sma- N m cosa) 


(A.9) 


(A.10) 


for the heavier Higgs, where a is the Higgs mixing angle. The quark coefficients are evaluated 
as 


with 



9 r (0 
2m w q 


r (i) 

' U 


for the up type quarks and 


r 


(i) 

d 


sin a 
shyd’ 


cos a 
sin ft 


cos a 
cos/3’ 



sin a 
cos ft 


(A.11) 


(A.12) 


(A. 13) 


for the down type quarks. The quark tensor contribution coupling in (2.15) can be expressed 
as 


<h 


E 


£ +b 2 ~ 

Qi Qi 


4 *=t 2 - {m^+m q ) 2 } 2 ' 

The gluon part of the scalar effective Lagrangian is fully determined by 


rriy rn~ 

b = —T~+ B d + B s - —B Zl 


25 


(B 1D + B 1S ), 


where 


i M 

t~ ~ 1 v z 9i 

9 96tt fy. 2 m 2 H . m 2 ~ 


Bn = 


E(«| - b l) m M m qn m q, m^) 


32vr 91 


Bo = 


32vr 


E(°| + b l m z ;) 


9.* 


^ = Toy E(«| - 


12vr 


9,* 


E(«! + &J)m Fi / 4 (m~, m^) 


9 ,* 


^25 = vx- E(«| + m^). 


487T 


9 ,* 


(A. 14) 

(A. 15) 

(A.16) 
(A. 17) 
(A. 18) 
(A. 19) 
(A. 20) 
(A.21) 
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The squark effective couplings in (B.16) are 


where 


c ( i } = 

91 


s (~^) (Tiq cos 2 9 q — e q sin 2 9 W cos 2 9 q 


cos 9 W 

.2 


C ( i } = 

92 


rriw 2mw 

s 0)(r 3 „ sin 2 9 q + e q sin 2 9 W cos 2 9 q ) 

cos 9 W 

+ - V, W ), 


m M / 


2mw 


s^ 1 * = — cos(a + /3)s ® = sin(a + (3) 


and 



cos a 

sin 13 ’ 
sin a 
cos/?’ 


„/( 2 ) 

u 

_/(2) _ 
r d “ 


sin a 
sin/3’ 
cos a 
cos (3 


(A.22) 

(A.23) 

(A.24) 

(A.25) 
(A.26) 


Loop integrals A - J 5 are given by Eqs. (Bla-e) in Ref. || (one must take care to correct 
the typo noted in Ref. f2|). 

For the only effective coupling needed in the spin dependent Lagrangian we have 


d q = 


V 


% + K 


4 . i „ m~ — (m 
*= 1 ,2 9i v 


z i 


+ m 9 ) 2 


4m^ 




3ij; 


where 

o"* = \(Nl - Nl) 

is determined by the neutralino mass matrix diagonalizing matrix N. 


(A.27) 


(A.28) 
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FIGURES 
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FIG. 1. A plot of contours of neutralino scattering events/kg/day in a ' 3 Ge detector, for 
mSUGRA parameters Ao = 0, tan/3 = 2 and a) fi < 0 and b) n > 0, in the mo vs. m i/ 2 
plane. The regions labelled TH are excluded by theoretical considerations, while the EX regions 
are excluded by collider searches for SUSY particles. The region to the right of the solid contour 
is excluded by tth 2 > 1. We also show contours of neutralino relic density f lh 2 = 0.15 and 0.4; the 
region in between is favored by models of a MDM universe. 
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FIG. 2. Same as Fig. 1, except for tan (3 = 10. Below the dotted contour is where flh 2 < 0.025 
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FIG. 3. Same as Fig. 1, except for tan/3 = 35. Note, however, the expanded scale relative to 
Fig. 1 and 2. 
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FIG. 4. A plot of neutralino scattering events/kg/day in a 73 Ge detector, for mSUGRA pa¬ 
rameters mo, mi/ 2 , Aq =150, 200, and 0 GeV versus tan/3 for /j, < 0 and /j, > 0. 





that from scalar interactions for rnSUGRA parameters mo,m 1 / 2 , Aq =150, 200, and 0 GeV, versus 
tan (3 for n < 0 and u > 0. 
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FIG. 6. A plot of neutralino scattering events/kg/day in a 73 Ge detector, for mSUGRA pa¬ 
rameters =150, 200 GeV with tan/3 = 10, versus the parameter Aq, for fi < 0 and 

u > 0 . 







FIG. 7. A plot of neutralino scattering events/kg/day in a 73 Ge detector for rnSUGRA param¬ 
eters Ao =150, 200, 0 GeV with tan/3 = 20 and fi > 0, versus month of the year, from 

Jan. 1 to Dec. 31, showing the expected seasonal variation in dark matter detection rate. 
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for LEP2 and Fermilab Tevatron MI experiments for detecting SUSY. 
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FIG. 9. Same as Fig. 2, but with in addition b —> 57 exclusion contours, plus the reach contours 
for LEP2 and Fermilab Tevatron MI experiments for detecting SUSY. 
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